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Abstract 



We analyze the steady-state flow as a function of the initial density for a class 
of deterministic cellular automata rules ("traffic rules") with periodic boundary 
conditions [H. Fuks and N. Boccara, Int. J. Mod. Phys. C 9, 1 (1998)]. We are 
able to predict from simple considerations the observed, unexpected cutoff of the 
1/-^ I average flow at unity. We also present an efficient algorithm for determining the 

' exact final flow from a given finite initial state. We analyze the behavior of this 

algorithm in the infinite limit to obtain for TZm,k an exact polynomial equation 



, maximally of 2(m + A;)th degree in the flow and density. 
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1 Introduction 

^ There is considerable interest in modeling traffic behavior via one- dimensional cellu- 

lar automata (CAs). The original models by Nagel and Schreckenberg[10] and Fukui 
and Ishibashi[6] are analyzed in e.g. [12, 13, 4], and the more general behavior of 
sum- conserving CAs is considered in [2, 1]. In [5], Fuks and Boccara introduced an 
interesting class of generalized deterministic traffic rules 7lm,k, which display a surpris- 
ing steady-state behavior: the expected flow of the cars never exceeds one regardless 
of the constraint values m and k. 

In these rules, as is usual for traffic rules, the road is represented as a one-dimensional 
lattice where each site has as its value either (empty) or 1 (car). Under 7lrn,k, a block 
of cars (ones) at most k units long moves right at most m units, or to the beginning 
of the next group. The same rule can also be expressed as follows: at each turn, each 
maximal match of 1^0^ is replaced (see Fig. 1): 
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road representation 



symbol representation groups 



01001111000000000010000110111101111000001 *_2,-i 1-1,2 07,-1 Oi,o 1-2,2 1-2,2 02,-1 7 

10011100011000000000010011110111100011000 l_i,i *o,o 07,-1 1-1,2 1-2,2 *o,o *o,-i 7 

00110001100011000000000111011110001100011 *-i,o *o,o *o,o 06,i 1-2,2 *o,o *o,o V 

11000110001100011000000101111000110001100 +0,0 *o,o *o,o 03,-1 1-2,2 *o,o *o,o *-i,o 8 

00011000110001100011000011100011000110011 *o,o *o,o ^0,0 ^0,0 Oi,i *o,o ^0,0 *-i,o 8 

01100011000110001100011010001100011001100 *o,o *o,o *o,o *o,o *-2,-i *o,o *o,o *-i,o *o,o 9 

Figure 1: An example time evolution under rule 7^3,2 from an initial state with L — 
41 into a cyclic state, expressed both in the road representation and the symbolic 
representation introduced in Sec. 4 in the text. Periodic boundary conditions are used. 

where a — min{A:,x} and h — min{m,t/}. From this representation, the dualism 
between the motion of the cars under the rule 'K-m.k and the motion of the empty sites 
under rule 'R,k,m in the opposite direction, as mentioned in [5], is obvious. 

The "physical" quantities of interest in systems that obey these rules are p, the 
density of ones, and 0, the flow, defined as piv)^ where (t)) is the average velocity 
of the cars. For finite-length systems, we write for the time-averaged steady-state 
fiow from a single state and for the average of over all states. For infinite-length 
systems, is the steady-state fiow. The equation 

expresses one consequence of the dualism discussed above. There are also other quan- 
tities such as acceleration, but these are outside the scope of this article. 

In this article, we examine the steady-state fiow of TZm,k: obtaining an exact polyno- 
mial equation in the infinite case. In the following sections, we first develop a formalism 
based on representing the road as a sequence of blocks rather than single sites. We 
show that the average fiow is fully determined by the number of these blocks in the 
steady state. In Sections 3 and 4, we use this fact to obtain simple upper and lower 
limits for the average fiow and present an efficient algorithm for calculating the steady- 
state fiow from a given finite initial state. In Section 5, we consider the behavior of the 
algorithm in the infinite limit and derive a steady-state condition, which we then solve 
in Section 6, yielding an analytical solution in the case of an infinite space. Finally, in 
Section 7, we obtain a non-trivial upper limit for the expected average fiow in a finite 
space. 

2 Fundamental properties of lZm,k 

The fiow of cars under rule TZm.k is easier to understand if the state of the road is 
considered as a sequence of blocks instead of single cars. As we shall see later, it is 
practical to distinguish between short, just, and long blocks, comparing the length of a 
block with m or /c as follows: a block of zeroes less than m sites long is a short block, 
more than m sites long is a long block, and exactly m sites long is a just block. For 
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Len; 


gths 


Po- > 0, po+ = 


Po- > 0, po+ > 


po- = 0, P0+ > 


Pi- 


> 0, = 


Cyclic: Intermediate 


Uncyclic 


Cyclic: Free-flowing 


Pi- 


> 0, px+ > 


Uncyclic 


Uncyclic pc may increase 


Uncyclic, pc may increase 


P\- 


= 0, pi+ > 


Cyclic: Congested 


Uncyclic, pa may increase 


Uncyclic, pa may increase 



Table 1: The different types of states in tlie system. As discussed in tfie text, tfie 
existence of short and long blocks distinguishes the different types of states. 



blocks of ones, the length is compared with A; in a similar fashion. Wc say that a pair 
of a 0-block and a 1-block is a group and define pc as the density of these groups. We 
also define po+ and po- as the densities of long and short 0-blocks, respectively, and 
similar symbols for the 1-blocks. 

The states of the system can be divided into nine different categories by the existence 
of short and long blocks, see Table 1. In the following, we will consider the three cyclic 
types of states separately. To verify that these are the only cyclic states we first show 
that the following two cases are unstable: long and short blocks of one kind, pi+ > 
and pi_ > 0, and long blocks of both kinds, pi+ > and po+ > 0. The other cases 
follow from the same proofs by duality. 

First, we note that new long blocks can never form, because a non-long block of 
cars moves continually to the right and therefore can not absorb other cars from the 
left. This also imphes that a non-long block can only be absorbed to an already long 
block on the right. The duality proves the same for long 0-blocks. Furthermore, the 
length of a long block can not grow, because long blocks emit just blocks whereas they 
can only absorb short or just blocks. 

To show that the states that have both short and long blocks of one type are 
unstable, consider the sequence 

where x < k, y > k, and 0-"* represents a 0-block that is either short or just. On 
each of the z first steps, the block 1^ absorbs one just block and emits one just block 
from the other end - its length remains unchanged. However, on the {z + l)th step it 
absorbs a block of length x and emits a block of length k. Therefore, the number of 
short blocks has decreased by one, and the length of the long block has decreased by 
k — x, possibly transforming it into a just or short block. Mathematical induction using 
this argument shows that, if all 0-blocks are short or just, then during the simulation, 
either the number of short or long 1-blocks drops to zero. The number of groups does 
not change in this process. An important observation is that the long block behaves 
like a decaying quasi-particle that is moving in the opposite direction from the ones 
by continuously absorbing short or just blocks and emitting just blocks. Naturally, 
applying the dualism property proves the same for long and short 0-blocks. 

Next, we show that if there are both long 1-bIocks and long 0-blocks, the state is 
unstable. We can first apply the above property to show that either a long 1-block 
decays or it eventually meets a long 0-block moving in the opposite direction (or else 
there are no long 0-blocks left in the system). But when the long blocks meet, they 
react and annihilate each other partially or wholly: the group 

on^ (2) 
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Description 




Conditions 









Free-flowing 


Pi+ 


= 0, 


Po- -- 


= 


mp 




Intermediate 


Pi+ 


= 0, 


Po+ -- 


= 


p{l- 


- p)Ipg 


Coiigostod 


Pi- 


= 0. 


Po- - 


= C) 




-p) 



Table 2: The three different types of cychc states for TZm,k and the formula for the flow 
(f) in each. 

where x > m and y > k emits just 0-blocks leftwards and just 1-blocks rightwards at 
each time step, reducing to 

and increasing the number of groups by one. Therefore, eventually either the long 0- 
blocks, the long 1-blocks or both will be exhausted, which shows that the initial state 
is not cyclic. 

This completes the study of the uncychc states, showing that from any initial state 
we will finally end up in one of the remaining states shown in Table 2. We discuss 
these three cyclic types of states separately below. 

It is fairly easy to sec that if there are only short and just blocks, then the 1-blocks 
move in the positive direction and the 0-blocks move in the negative direction but the 
number of blocks and the distribution and relative order of the and 1-blocks among 
themselves do not change - the state is obviously cychc. To evaluate the fiow, we first 
note that in such a state each block of cars travels on each step on average (1 — p)/ pa 
units forwards, which, when multiplied by the density of cars yields 

= P0+ = 0, pi+ = 0. (4) 

PG 

As this class of states does not correspond to either the free-fiowing or congested 
phases of the simpler traffic rules TZ^^k and 7?.m,i, we term it, for want of a better name, 
intermediate. 

The free-flowing states where all 0-blocks are long or just and all 1-blocks are short 
or just are also simple. All the cars obviously move forwards at maximum speed and 
consequently, these states are also cyclic, with 

= mp, po- = 0, pi+ = 0. (5) 

Applying the dualism between zeroes and ones we obtain the formula 

= A;(l-p), po+ = 0, pi_ = 0, (6) 

for the opposite case: congested states with long or just 1-blocks and short or just 
0-blocks. 

We can summarize the above by noting that p and the final pa determine the type 
of the cyclic state. This follows from the fact that the final average lengths of 0- 
and 1-blocks can hold for only one type of a cyclic state: in the intermediate phase 
blocks must on the average be short or just whereas in the free-flowing and congested 
phases either 0- or 1-blocks must be long and the other blocks short or just as shown 
in Table 2. Writing these conditions in terms of p and pa allows us to combine the 
above evaluations of (j) into one surprisingly simple formula: 
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Proposition 1 The average flow over a cycle in any cyclic state of 7lm,k is 

(f) = min{mp, p(l - p)/pG, k(l - p)}. 

The treatment of long blocks above also yields the following proposition. 

Proposition 2 During evolution of the system, pa can never decrease. It can increase 
only when po+ > and pi^ > 0. 

These two propositions give the system its interesting characteristics: the final flow 
is completely determined by the final value of pa, which in turn depends on the intricate 
reactions of the long blocks. 

As a direct consequence of these propositions, it is straightforward to obtain crude 
upper and lower limits under both finite and infinite length (see Fig. 2): 

Proposition 3 The flow of any cyclic state ( and thus the average flow over different 
states) satisfies 

> min{mp, |p- 1/2| + l/2,/c(l -p)} (7) 
(f) < minjmp, A;(l — p)} (8) 

for any < p < 1 and any size of the lattice. 

The lower limit is obtained through considering what is the greatest possible number 
of groups. For a lattice of size L, this is easily seen to be 

which together with Prop. 1 yields the first part. The second inequality follows directly 
from Prop. 1. 

When m — 1 or k — 1, this proposition reduces to the well-known flow formula 

= min{mp, k{l — p)} 

discussed e.g. in [5]. 



3 Simple upper and lower limits at infinite length 

At inflnite length, we can use block probabilities 

P(aia2 ■ • • a„) = pS.K=il(l - p)E.K=o] (9) 

in the initial random state in order to calculate various statistics. The brackets [■ ■ ■] 
represent Ivcrson's notation[8, p. 11], which evaluates to 1 if the enclosed statement is 
true and to otherwise. 

The frequency of groups in the initial, random state is obviously given by the density 
of group edges: 

PG,imtmi=P(Ol) =P(1-P)■ 



(a) 



trivial upper 
trivial lower 
upper 
lower 
exact 




(b) 




Figure 2: Different limits and the exact solution for the flow as a function of density 
p at infinite length for (a) 7^2,2 and (b) 7^3,2- The trivial limits are given by Eqs. (7) 
and (8). The equations (11) and (12) yield the tighter limits. The enclosed regions are 
shown magnified in Fig. 7. 



Proposition 2 tells us that 

PG > Pg, initial- 

Combining these with Proposition 1 yields the important upper limit 

< min {mp, 1, — p)} , (10) 

which is valid for any TZm.k, as Fuks and Boccara observed experimentally in [5]. The 
expected flow is cut off by the fact that initial groups can never combine into fewer 
groups, which limits the maximum speed of cars. We shall show in Section 7 that this 
limit is also valid for finite-length systems. 

This hmit does not take into account the dynamics of the system but assumes 
that the final state has the same number of groups as the initial state. It is possible 
to improve this upper limit by explicitly including some cases where certain initial 
configurations are known to produce new groups. For adjacent long blocks of zeroes 
and ones that directly react with each other as shown by Eqs. (2) and (3), we have 



PG > E HiinjW Jfj}p(10n^0) 

x,y>l ^ ^ ^ 



r>0 l-p'=(l-p)-' 

which yields the tighter upper limit 

< min{mp,l -^^=(1 -p)"',A;(l -p)} (11) 

when combined with Proposition 1. 

A lower limit can be estimated by computing the density of groups that could arise 
by having all long blocks split maximally. For ones, this is 



X>1 



r>0 ^ P 

Combining with the same formula for zeroes yields 

(j) > minjmp, max{l - p^ 1 - (1 - p)™}, k{l - p)|. (12) 

Figure 2 shows the accuracy of these limits for TZ2,2 and 7^.3^2 ■ For larger m and 
k, the limits become tighter, converging to unity geometrically as m and k tend to 
infinity. 



4 Efficient algorithm for computing 
the steady-state flow 

It is not necessary to carry out the simulation to find the number (or density) of 
groups in the final state, and thus the steady-state flow. In this section, we present 
an 0{L) algorithm for finding the final number of groups from a given starting state 
with periodic boundary conditions in a lattice of length L. This is interesting from 
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two different perspectives: first, it makes it possible to calculate the final ffow for long 
strings more effectively. Additionally, it helps us to understand the dynamics of the 
system and derive analytic results on the behavior of the system in the next Sections. 

Let us define the symbols -ka,b, '^a,b, ^a,b and la,b all to correspond to 0"+*^!^+"* for 
different a and h: 

*a,fe, a < and 6 < 0, 
b+m ^ I 0a,6, a > and 6 < 0, 

la,6, a < and 6 > 0, ^ ' 

'^a,bi a > and 6 > 0. 

We represent the initial state in this notation and then carry out a series of string 
replacements. The final number of groups is obtained as the initial number of symbols 
plus the number of extra groups created by the replacements as shown below. 
Diamonds react as follows: 



Oa,b ^'^a-m,b-k, ^9 = +1, (14) 

where the wildcard symbol ? represents the correct symbol for the new indices from 
Eq. (13) and /S.g represents the number of new groups created by the reaction. 

Stars can be collected, along with zeroes and ones, but zeroes only from the right 
and ones only from the left. No new groups are created. 

= 0, (15) 
= 0, (16) 
A^ = 0. (17) 

Finally, when a and a 1 meet, it is possible to form new groups: 

^aAc,d ^?a+c,6+d, = 0. (18) 

Note that new groups are not created directly by this reaction, but the result may be 
a O. 

The above replacements recursively evaluate the reactions between long blocks. To 
prove the correctness of a maximal application of the replacements, we consider the 
reactions occurring in an actual simulation. As described in Section 2, all new groups 
arise from the annihilations of long blocks and the length of a long block can never 
grow. Thus, the number of extra groups is dictated by the sequence of short blocks that 
reacting long blocks must absorb before annihilating. In the following, we shall show 
that each of the nested annihilations in an actual simulation is correctly represented 
by the string replacements. 

First, consider a subsequence starting with a long 0-block and ending with a long 
1-block with only short and just blocks between them. Suppose further that the long 
blocks are long enough to actually meet before turning into short or just blocks. Then 
a maximal application of Eqs. (15)-(18) to the subsequence subtracts the total "short- 
ness" of the short blocks from the long blocks yielding a diamond symbol, which cor- 
rectly represents what remains of the long blocks when they meet. The diamond will 
then annihilate according to Eq. (14) and yield the correct number of extra groups 
and a symbol representing the remainder (cither a ^, if both long blocks are completely 
annihilated, or a or 1, if the annihilation is partial). This remainder symbol is exactly 



^a,b'^c,d ^ ^a+c,b+d: 
^a,b~^c,d ^ ?a+c,6+<i) 
'^a,b^c,d ^ • a+c,b+d: 



what would be the result of an actual simulation of the subsequence minus the just 
blocks that the sequence would have emitted from both ends. Because the just blocks 
do not affect long blocks outside the subsequence, they can be disregarded. 

Suppose then that either long block would decay before meeting the other long 
block. Then an application of Eqs. (16)-(18) to the subsequence will at some point 
turn either the 0- or 1-symbol into a -k (or the whole sequence can turn into a 0, 1, or 

A long block at one or both ends of the subsequence has thus turned into a just 
or short block, making it possible for outside long blocks to react over the remaining 
just and short blocks of the subsequence. Again, the just blocks that a decaying long 
block would emit are disregarded as they do not affect the length or order of other long 
blocks. 

These are in fact all the cases that need to be considered, as the above cases can 
be applied recursively to the results of inner annihilations. If neither case applies, we 
have found the final number of groups, because either long 0-blocks or long 1-blocks 
have then been exhausted and so there cannot be any other reactions in an simulation 
nor can the replacements form any new O's. On the other hand, each of the applicable 
reactions will be carried out at some point of a maximal application of the replacements, 
because the and 1-symbols at the ends or a O do not interact with outside symbols. 
The replacements may evaluate the decay of the annihilating long blocks in different 
order from the actual simulation, but as the long blocks must absorb all short blocks 
between them before they can meet and the short blocks can not disappear unless 
absorbed by a long block, the result is still the same. 

Note that the remainder of a O may be any non-O symbol, causing complicated 
recursion as the result of an annihilation affects further replacements. For example, 
consider the replacements 

1 * 



where the innermost reactions must be evaluated before the outer reactions can be 
considered regardless of the order of the possible replacements. 

The iteration of Eqs. (14)-(18) can be carried out by the stack algorithm shown in 
Fig. 3. Each symbol of the input string requires 0{l) operations to be dealt with in this 
algorithm, where / is the bit-length of the symbol. This also includes the final steps to 
wrap the stack. Therefore, the running time of this algorithm is obviously 0{L) for a 
lattice of length L. 

The worst-case time for simply running the cellular automaton simulation, on the 
other hand, is ^(L^) since in this system far-away cars do interact with each other. 
For example, the initial state 0'"+^(1^0™)^^^l^+^ (in the road representation) would 
require T simulation steps to reach a cyclic state. 

Figure 4 illustrates an apphcation of the algorithm on an example string. 



0. Initialize and create an empty stack. 

1 . If there are no symbols left in the input string, go to step 6. 

2. (Input new symbol) Read the next symbol and push it to the stack, and setg ^ g+1 

3. (Reaction to create new groups) While the top symbol is a O, i.e. (a, b) with a > 
and b > 0, subtract (m, k) from the symbol, and set ^ + 1. 

4. (Drop reaction result down stack) If the top symbol is a non-0 and the symbol 
below it is a 0, pop the two symbols off the stack and push their sum and go to 
step 3. 

5. Go back to step 1. 

6. (Finish up by wrapping the stack) If the top symbol is a non-0 and the bottom 
symbol of the stack is a 0, remove the bottom symbol and push it back to the top 
of the stack, and continue from step 3. 

7. Terminate. The variable g now contains the total number of groups that this string 
will have in the steady state. 



Figure 3: This algorithm calculates the final number of groups in the given symbol 
string with periodic boundary conditions. The symbols 7a,b are represented as pairs 
(a, b) and the input string is read left to right. Addition between pairs is defined as 
(a, b) + (c, d) = {a + c,b + d). Step 3 corresponds to Eq. (14) and steps 4 and 6 both 
correspond to Eqs. (16) and (18) in the text. 



in 



time groups dropped symbols active stack remaining input 








*-2,-l 1-1,2 07,_i Oi,o 1-2,2 1-2,2 02,-1 


1 1 




^—9 —1 




1-1,2 07,-1 Oi,o 1-2,2 1-2,2 02,-1 


2 2 -*--2.-i 




— 1 ^'1 




O7 -1 0-| 1-0 ■-> 1-0 ■-> Oo -1 


3 3 2 —1 1—1 2 




f \ 
07,-1 




Ol,0 1—2,2 1—2,2 02,-1 


4 4 *_2,-l 1-1,2 




0? 1 

^7,-1 


Oi,o 




1-2,2 1-2,2 02,-1 


5 5 •*--2,-l 1-1,2 
5 *-2,-l 1-1,2 

5 *-2,-l 1-1,2 

6 *-2,-l 1-1,2 




07,-1 
07,-1 
06,1 
03,-1 


Oi,o 


,1-2,2 J 


1-2,2 02,-1 
1-2,2 02,-1 
1-2,2 02,-1 
1-2,2 02,-1 


6 7 *_2,-l 1-1,2 

7 *-2,-l 1-1,2 

8 *-2,-l 1-1,2 




03,-1 

Oi,i 

*-2 -1 


1—9 9 




02,-1 
02,-1 
02,-1 


7 9 *-2,-l 1-1,2 *-2,-l 




^ \ 

02,-1 






9 1-1,2 *-2,-l 
9 1-1,2 ^-2,-1 




02,-1 
.*0 -2 . 


■>^-1 -1 







Figure 4: Evolution of the stack of the algorithm in Fig. 3 operating on the sample 
string in Fig. 1. In the end, the number of extra groups is the same as in Fig. 1. The 
doubled line signifies the end of the input string, after which the stack is wrapped over 
(cf. step 6 of the algorithm in Fig. 3). The division to dropped symbols and active 
stack is implicit ~ it is only significant for the Markov model analysis of the algorithm 
in Section 5. The boxes depict the "lifetimes" of the sub-stacks started by each input 
symbol. The time steps of the Markov model are at the bottom fines of these boxes, 
where the sub-stacks are finished. 
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5 Steady-state conditions for the stack machine at 
the infinite hmit 



It is often the case that infinite hmits of systems are easier to solve than finite cases. 
This is also the case with TZm,k' in this section, we examine the behavior of the stack 
machine algorithm as the length approaches infinity. We regard the evolution of the 
stack configurations as a Markov process and derive a set of equations for the sta- 
tionarity of the probabilities of stack configurations. These equations are solved in 
the next Section to determine the probabilities of different annihilations and thus the 
steady-state fiow. 

Consider the symbols and reactions defined above in Section 4. All new groups are 
created from the diamonds, which in turn can only arise when a 0-symbol is combined 
with a following 1-symbol. The algorithm in Fig. 4 uses this fact to find the final 
number of groups by only tracking the reactions of 0-symbols. It scans through the 
input string linearly, from left to right, maintaining a stack of the processed symbols 
with all reactions of O's and O's carried out. This means that all the remaining 0- 
symbols end up on top, which we shall now call the active part of the stack. When a 
non-0 input symbol consumes all the O's, we say that the symbol is dropped off from 
the bottom of the active stack into the inactive part, as the symbol can no longer create 
new groups with the following symbols. 

In finite systems with periodic boundary conditions, the processing of the input 
string is divided into two parts. First, the input string is scanned as above. When 
all of the input has been read, the inactive part of the stack, comprising of I's and 
T*r's, is reprocessed with the zeroes in the active part, since the I's on the bottom of 
the stack may react with the O's on top, producing new groups. The relative effect of 
this wrapping diminishes as the length approaches infinity: it is easy show that we can 
ignore all symbols dropped to the inactive part in the infinite limit. 

To obtain the limiting flow, we thus need to evaluate the average number of new 
groups produced as the stack algorithm processes a new symbol. When new symbols 
are input, the active stack can either grow infinitely or remain finite. If the active 
stack remains finite, the probabilities of different active stacks will eventually reach a 
stationary distribution. Once the stationary distribution is known, it is straightforward 
to calculate the expected number of new groups for a random input symbol. 

If, on the other hand, there are not enough I's and -k^s to annihilate the O's and 
the stack grows infinitely, we can use the dualism property and consider the thus finite 
stack of I's instead of O's. It turns out later that we do not even have to consider this 
dualism explicitly, as the symmetry of the equations is restored in the next section. 

Formally, we regard the evolution of the active stack as a Markov process. The 
Markov property is satisfied as the next state depends only on the current state and 
the upcoming independently distributed random symbol. Clearly the process can reach 
each possible stack configuration from every state and has a positive probability to 
remain in its current state (the symbol -kQ ^). Furthermore, given that the active stack 
will not grow infinitely, the process will return to every state an infinite number of times 
with probability one. This means that the Markov chain is irreducible, aperiodic, and 
Harris recurrent and therefore will converge to a unique stationary distribution (see, 
e.g., [3] or [11, Proposition 6.3]). 
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An essential property is that the algorithm can be applied independent of the lower 
levels for each stacked symbol until that sub-stack is finished, that is, until the lowest 
level of the sub-stack turns into a 1 or a which happens either immediately, if the 
stacked symbol is already a 1 or a or when a 1 or ★ on higher level falls to the bottom 
level consuming all O's on the sub-stack. This allows us to consider each level of the 
active stack as the bottom of an identically distributed sub-stack. The distribution 
is particularly interesting when a sub-stack is just finished. At these times the active 
stack consists of zeroes at the bottom of each sub-stack and of the 1 or ★ that finishes 
the topmost sub-stack (see Fig. 4). 

Suppose p{l.x.y) is the distribution of symbols seen on the bottom level of the ac- 
tive stack at each time a sub-stack is finished. Then, if the symbol is a 0, the same 
distribution of symbols will be seen on the bottom of the sub-stack above that symbol. 
Thus, a symbol distribution defines a distribution for stack configurations. Note that 
an arbitrary stack distribution can not be represented by such a symbol distribution 
but it is required that the stack symbols arc identically and independently distributed 
and that the height of the stack is implied by the stack symbols as described above. 
Furthermore, even though each input symbol starts a sub-stack so that there is the 
same number of time steps as there are input symbols, the sub-stacks are finished 
out-of-synch with the time steps of the algorithm. This complication, however, is in- 
consequential as the expected density of new groups on each time step is still the same 
as the expected per-symbol density. 

For simplicity, we consider the input symbol distribution Pi{'?x,y) to represent what 
remains of the symbol after initial reactions of O-symbols. This difference is only 
conceptual as the algorithm would immediately carry out the initial reactions for each 
input symbol anyway. It is easy to see that the resulting distribution must still be 
geometric; only a constant is required to normalize the lack of O's. The normalized 
input distribution is defined for indices in the set 

D -.^ {{x,y) : X > -m, y > -k, [x > 0][y > 0] = 0} 



as 

where we define the symbol A to represent the quantity 

A := (l-p)V, (20) 

which occurs often in the formulas below. This quantity is the initial probability of a 
O, and the initial reactions produce a density of A/{1 — A) new groups (cf. Sec. 3). 
Here and in the following, we use densities relative to the initial symbol density 1 in 
the string of symbols. 

In our model the initial stack distribution is defined by p{-\t = 0) := Pi{-) cor- 
responding to the distribution of stacks that results from running the algorithm on 
random input until the first non-0 symbol is stacked in, i.e., at the first time a sub- 
stack is finished (cf. Fig. 4). 

The transition from a symbol distribution p{-\t) to the distribution p[-\t + 1) on 
next time step can be defined by considering the possible ways for a given symbol to 



arise on each level of the stack: A symbol can be the result of a and a non-0 above 
it reacting as per Eqs. (16) and (18). The reaction will result in a non-O symbol lx,y 
with probability 

Pc(?x,3/|^) = Y.Pi^^-'^A't)p{'^-a,y-b\t). 

a<x; a<0 b<0 

If, on the other hand, the result of the reaction is a diamond, it will further react, 
possibly several times, according to Eq. (14) and yield 'Ix^y, where {x,y) G D, with 
probability 

Pd{^x,y\t) = YYYp{Oa:-a+rm,b\t)p{la,y-b+rk\t). 
r>l a<Ob<Q 

In case there is another on top of a Oa,^^, no reactions will occur and the O^^y remains 
for the next time step. Finally, if the symbol is not a 0, it falls off from the bottom of 
the sub-stack and is replaced by a fresh symbol from the input distribution. Thus, the 
transition function can be written as 

= [{x,y) e D]pd{l^Jt) + [ix,y) G E]p,{?,Jt) 

+ [x> 0]p{7,Jt)pm + (1 - p(0|t))pi(?,,,), (21) 

where p{0\t) denotes the total probability of 0-symbols and the set E is defined as the 
complement of diamonds: 

E:^{{x,y):[x>0][y>0]^Q}. 

Note that the transition function does not define a Markov chain for symbols but it 
implicitly defines a linear Markov operator for the subspace of stack distributions that 
are determined by symbol distributions. 

The stack distribution defined by a symbol distribution is stationary when 

P{'^.,y\t + l)^p{l,,y\t). (22) 

Thus, we need to find a symbol distribution p{'fx,y) corresponding to the unique sta- 
tionary stack distribution. The symbol distribution can then be used to determine 
probabilities for different reactions. 

Figure 5 represents the possible indices for different stack symbols. It is easy to see 
why some symbols can not occur at all, but even more can be said. The distribution 
of the stack symbols retains some of the geometric properties of the input distribution. 
The distribution of 0-symbols is geometric in its x index and the distribution of 1- 
symbols is geometric in its y index. (The distribution of Y^. is in fact geometric in y 
for all {x,y) G D.) This can be justified by noting that two such symbols result from 
the same set of paths of the process with the corresponding difference only in a; or |/ 
index of one specific input symbol (for O's, the initial starting the sub-stack and for 
I's, the final 1 that falls through to the bottom of the sub-stack). These properties are 
listed below: 

Pi'^x,y) = 0, X < -m, 

p{Ox,y) = 0, x>0,y>0, , . 

p{0,,y) = (l-p)=^-V(OiJ, x>0,y<0, ^ ^ 

p{ix,y) = pf-V(l,,,i), x<0,y>0. 



14 



D 





— m 


short just 
1 — m 2 — m ■ ■ ■ 


long 

1 2 


-1-k 
short —k 
1-k 

just 


•-. ^ \ N \ 
{*-m.,-l-k) 


\ \ \ 


\\\\\\-\\ 
) *2-m,-l-fc • • • *0,-l-A; 


\ ^ ' 

Ol,-l-A; 


^ \ \ \ \^ 

02,-1-k -x 












1 

long 2 








(^m) (^2,i) ■■■ 
(Ol,2) (02,2) ■ ■ ■ 



Figure 5: The possible index combinations for different symbols in the Markov model 
for the algorithm discussed in Sec. 5. By considering the possible outcomes of combined 
and non-0 symbols, it can be seen that the parenthesized symbols can not occur. The 
sets D and E defined in Sec. 5 are depicted by the different shadings. The marked 
sequences are represented by the generating functions Fi_rn{z) and Fi{z) discussed in 
Sec. 6. 



For a more rigorous proof, it is easy to check that the initial distribution pj(-) has 
these properties and that the transition function maintains the properties. Thus, the 
limiting stationary distribution must also have the properties. 

The stationarity recurrence given by Eqs. (21) and (22) can be transformed to 

{l-[x> Ob(0))p(?,,,) _ [{x,y) e D]p,{l,,y) 



(1 -p{o)){i-p)-py {i-p{o)){i-prpy 

■x,y) Pi[-x,y) /r,A\ 

(i-p(o))(i-p)v^ {i-p^py' ^ ^ 

For clarity, we define f^^y as the left side of the equation: 

:= (l-P)-V'i^J^|«p(?.,). (25) 

This transformation of p{-) cancels out the geometric factors and will decouple the re- 
cursive p{0) coefficient from the stationarity equation. The transformation is reversible 
as p(0) can be obtained in terms of fx,y from 

p(0) = ^ p{Ox,y)^ Y.i^-Prp'Uy (26) 

x>0,y x>0,y 

With this definition, the geometric properties reduce to 

fx,y = fx,i, x<0,y>0, ^27) 
fx,y = fi,y, x>0,y<0. 



We define analogously the components d^^y and c^^y corresponding to terms on the right 
side of Eq. (24) and apply the definition of fx,y and the above properties: 

V y l-p(O) 

— ^ ] ^ ^ ! ^ ! fx—a+rm,b fa,y—b+rk 
r->l^ o<0 6<0 

~ 1 _ A X/ /l''' 
^ a<0 6<0 

for (x, y) E D and 

^ ^ ^ ^ fx—a,b fa,y—b 

a<0 6<0 

— X^ X^ /l,6 fa,y-b 
a<Ob<0 

for (x, I/) e ii^. The stationarity condition given in Eq. (24) can then be expanded as 

fx,y = [{x, y) e D]dx^y + [{x, y) e E]cx^y 
+ (l-p)-V>i(?.J 

^ a b 

+[{x,y)eE] Yl Y.h,bfa,y-b 

a<x\ a<0 b 

where we have left out the summation limits for zero terms, based on f^fi = for a > 
and 6 > 0. 

Noting that only the middle line really depends on x and y and that fa,- — for 
a < —m, the recurrence can be written as 

fx,y = [{x,y)eD]- 

+ X] X] /l>ft fa,y-b 
a<x,a<0 b 

- [x>0][y>0](^-l), (28) 



C 



where we define 



/l-m,l 2^ 2^ h,b Ja,l + 1 
a 6 

for reasons to become clear later. 

We have thus reduced the stationarity of the stack distribution given by Eqs. (21) 
and (22) to the above convolution equation, where fx,y determines p{'^-x,y) and is given 
by Eq. (25). 



Ifi 



6 Exact solution for the steady-state flow through 
generating functions 

The convolution recurrence in Eq. (28), which is the stationarity condition for the 
stack distributions, can be solved using generating functions (cf. [7]). We define a 
formal generating function Fx{z) for the sequence fx,, by 

F.i^) := EUyZ'. (30) 
y 

This generating function is not quite ordinary: the sum goes over all y, positive and 
negative. In general the values of a generating function do not uniquely determine a 
sequence that is positive at an infinite distance in both directions. Here, however, we 
know that both Fi{z) and Fi^rn{z) are uniquely defined generating functions, because 
when X = 1 — m, the term f^^y is only positive at an infinite distance in the positive 
y direction and when x = 1, it is only positive at an infinite distance in the negative 
y-direction (see Fig. 5). The coefficients for the functions for other x are positive 
infinitely in both the positive and negative y-dircction but they arc only used in the 
following calculations for formal multiplication and addition operations corresponding 
to well-defined convolution and sum operations on sequences. 

The term /a, ^ given by Eq. (28) vanishes for x < —m. Thus, we can represent it 
by the first non-zero case /i-m.y = — k]A/C, and the differences 

fx,y ~ fx-l,y — ^ 1] fl,b fx-l,y-b 

b 

-[x = l][y>Q]^^ (31) 

ioT X > 1—m. With the generating function notation we have Fi_^ (2;) = {A/C)z^~'^ /{l- 
z) and 

Fx{z) - Fx_,{z) ^ F,{z)Fx^,{z) ■ 1-^[^ = 1> 



for 1 — m < a; < 1. Thus, 



C l-z 



^x+m-l TT, /.^ 1 — C [x — l]z 



Fx{z) = {F,{z) + ir--'F,_^{z)- ^ (32) 

for — m < X <1. Thus, if we can solve Fi{z) and C, we have determined Fx{z) for all 
X, because C determines Fi^„i{,z). 

For X = 1, equation (32) can be written as 

(F,w + i) = (f,(.) + ir^^-^^, (33) 

which is an mth degree equation with respect to [3 := Fi{z) + 1. It is easy to see that 
there are at most two positive real solutions for /5. 

Because we can solve Fi[z) given C, the complete solution for the stationary dis- 
tribution of the stack configurations and thereby <p now hinges on determining C. 
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Unfortunately, the quantity C cannot be solved directly from the above equations, 
since its definition is already used in solving them; equations relating C and Fi{z) 
reduce to identities when combined with Eq. (33). 

However, there is a different, strange approach: we can determine C based on the 
fact that Fi(z) represents (indirectly) a probability distribution. The correct Fi{z) 
must obviously be analytic for some region \z\ > r. Additionally, it must be positive 
and decreasing in z, because it has non-zero coefficients only for non-positive exponents 
of z, and all coefficients must be nonnegative since they are probabilities multiplied by 
a positive function of the index. These two constraints allow us to uniquely determine 
C in the following. 

Eq. (33) can be solved with respect to C as 

(note that (3 = Fi{z) + 1 depends on z). Substituting a := in this equation yields 
a perfectly symmetric form 

^ - a(3-a-p- ^^^^ 

Figure 6 depicts the solutions of Eq. (35) with different values of C for 7^2,2 and 7l3,2- 
The figures are essentially similar for larger m and k with at most two positive solutions 
and in addition one everywhere negative solution if m is odd. In either case, it can 
be seen that a too large value of C results to a gap in the solution and a too small 
value results to either a non-monotonous or a non-positive solution. Only the correct 
C allows changing branches in the singularity point to obtain a feasible solution. This 
is analogous to the singular behavior of elliptic curves (cf. [9]). 

Because the surface C{a,(3) is smooth, the two constant-C contours can only cross 
at a critical point. The critical points 'VC{a,/3) = are determined by the equations 

f Aka''-^(3"'{a(3 - a - (3) - (Aa^/5" - l)(/3 - 1) = 0, 
\ Ama''l3'^-\ap - a - p) - {Aa'^p'^ - l){a - 1) = 0. 

Multiplying by {a — 1) and —{(3 — 1) and adding, the equations yield k{a — l)/3 = 
m{P — l)a. Changing variables to a = (a — l)/{ma) and b = {P — l)/{k(3) yields a 
simple form a = 6 for the equation. Substituting a = 1/(1 — am) and (3 = 1/(1 — ak) 
to the critical point equations and to Eq. (35) results in 

' A{l-am)-''{l-aky'^ 



l-a{k + m-l)' 

m 1 ^ l-a{k + m) Wt^J 

A(l — am] (1 — ak) =1 — C-, -i-, 

^ ' ^ > {I - ak){l - am) 

Now, we only need to reduce a from this system and then we have an equation relating 
the unknown C to the parameters A, m, and k. By eliminating the left sides we obtain 
from the right sides either a — l/{k + m) corresponding to a pole of Eq. (35) or 

C[l-a{k + m-l)] = {l-ak){l-am). (37) 
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Figure 6: Plots of the real solutions f3 = Fi{l/a) + l of the generating function equation 
(35) at p = 1/2 for C below, at, and above the correct (p for (a) 7^2,2 and (b) 7?.3,2. 
Because /3 is an ordinary generating function for a sequence of positive values, the 
correct solution and all its derivatives must be positive. Furthermore, because a — 1/p 
corresponds to a sum of probabilities, /? must converge at least for \a\ < 1/p. This 
excludes too large values of C. For too small values, the lower curve continues below 1 
yielding negative probabilities and the upper curve is not feasible either, because it is 
decreasing at a = Only the singular curve at the correct C yields an admissible 
solution. Because Eq. (35) is symmetric with respect to {a, k) and (/3, m), the solutions 
have symmetry axis a — (5 when m — k. 
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This is a second degree equation for a and its smaller solution is 



l + (l-C)(A; + m-l)-c 
^ 2k^ ' ^^^^ 



where we define 



c := ^J{l + {l- C){k + m- l)f - 4(1 - C)km. (39) 

Note that the other solution with +c in Eq. (38) does not yield a critical point. 
Now Eq. (37) can be used to rewrite the first equation of Eq. (36) as 

A = Ca{l-am)''-^{l-akr-^ (40) 
= C''a{l-a{k + m-l))''-^{l-akr-\ (41) 

yielding A as a function of C in closed form with the solution of a given by Eq. (38). 
In the special case of m = A;, we can solve p from A in closed form and obtain the 
density as a function of flow. 

For any given m and k, the equations (38)-(41) for A and C can be expanded 
to a polynomial equation, which is easily seen to be second degree in A and at most 
2(m + A;)th degree in C. In practice the degree of C seems to reduce to m + A; + 1, and 
we conjecture that this holds for all m and k. For example, in the case oi m — k — 2, 
the equation can be written out as 

16^2 + 8AC^ - 36AC^ + (1 + 27A)C^ - = 0. 

Now that we have determined all variables, we can determine probabilities for dif- 
ferent reactions. The density of new groups on the bottom level of the active stack 
is 

PAg = J2^J2J2p(^^-'i+rm,b)p{'i-a,y-b+rk) 

ix,y)eDr>l a<Oh<0 

= (1-P(0))E E^^''E(l-p)V/i,;./a,i 

(x,y)eDr>l a,b<0 

A 1 - A n - A \ 
= (1 -fC)) (137)5— (—-V- 

When all levels of the stack are taken into account, the total density of new groups is 

PA,(l+p(0)+p(0)^ + ...) = -^ ^ 



C 1-A 

Adding in the initial density 1 and the density of groups arising from the initial reactions 
yields the final group density relative to the initial group density: 



p(l-p) l-A \C I- A) C 

Thus, in the intermediate phase, (p — C. Furthermore, the phase transitions occur 
where C as a function of p crosses mp or /c(l — p), the fiow of free- flowing and congested 
phases. For example, when m — k — 2 the phase transitions can be solved to be exactly 
at 



2 7 V 2. 
The most important results are summarized below: 
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Proposition 4 The flow at infinite time for random infinite strings in the intermediate 
phase is C , where 1/C is the expected number of extra groups arising per each input 
symbol and is determined by Eqs. (38)-(41)- The flow depends on p only through the 
quantity A — (1 — pY^p^ . 

Proposition 5 For any given m and k, the flow (j) at infinite time and the density p 
in the intermediate phase can be related by a polynomial equation maximally of degree 
2{m + k). 



7 Upper limit for steady-state flow in flnite systems 

Carrying out calculations for finite systems is considerably more difficult, since the 
probabilities are no longer independent of each other. However, the following interest- 
ing limit can be derived. 

Proposition 6 The average flow in steady states starting from random binary strings 
of length L with < p < 1 satisfies 

(f) < min ^mp, ^ ~ ' ^(^ ~ P) 

where equality applies at least when k > L — 1 and m > L — 1. 

Proof. If < p < 1, the number of different initial states with a given number of 
groups can be counted by considering different ways of placing the group boundaries 
on the string. The distribution of Nq for random binary strings can be simplified to 



p{Ng) = < 



'pL-l\({l-p)L-l' 



Ng 



-,iiNG>0, 



0, if Ng < 0. 



Using basic binomial coefficient sum formulas (see e.g. [7]) and noting that Ng cannot 
be zero, the expected value of L/Ng is 



L 



Prom this, the formula in the Proposition follows. Finally, the equality follows simply 
from the fact that if the condition given holds, no groups can split. 

Note that with reasonable L and pL, the reciprocal of the binomial coefficient is 
negligibly small compared to 1. 
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Figure 7: Infinite-space behavior of fiow as a function of density under (a) 7^2,2 and 
(b) 7^3,2- Average of 100 simulations with L = 10000, theoretical flow, and upper and 
lower limits given by Eqs. (11) and (12). 

8 Simulations 

Simulations were carried out to test the theoretical results. For small L, complete 
summations were possible so the simulated curves are in fact exact. For large L, a 
number of random initial configurations were generated and the evolution of the system 
simulated. Since the resulting steady-state flow under TZm,k, when m > 1 and A; > 1, 
depends on the whole initial state (and not just p, as for when m = 1 or k = 1), the 
samples so simulated will in general not have the same flow. Therefore, the standard 
deviation is displayed along with the average of the resulting flows, giving an idea of 
the strength of the fluctuations. As L tends to inflnity, the fluctuations slowly average 
out, displaying the usual 1/VL behavior for the standard deviation. 

Figure 7 depicts the simulated flow and exact solution for infinite space. The 
theoretical solution agrees well with simulated results. Figure 8 shows the simulated 
fiow and upper and lower limits for finite space. For m, k > L — 1 the upper limit is 
exact as confirmed by the simulation. 

9 Conclusions 

In this article, we have solved the behavior of the generalized traffic rules TZm,k for 
infinite lengths of road and uniform random initial conditions. 

We have derived an efficient algorithm for computing the average fiow from an 
initial state under the generalized traffic rules TZrn,k- The idea behind the algorithm is 
an appropriate representation of the road as a string of 0^1^ blocks instead of single 
sites, and the fact that finding the average flow can be reduced to flnding the number 
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Figure 8: Finite-space behavior of the average flow as a function of density under 
(a) 7?.2,2 and (b) 7^7.7 when L = 8. AW points represent either evaluated formulas or 
averages calculated over all configurations. The trivial lower limit is given by Eq. (7). 
The upper limit given by Prop. 6 is exact for m, A; > L — 1 as demonstrated by b). 
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of these blocks in the cychc state. 

The algorithm works by decoupling the time from the simulation and considering 
directly the different reactions that would happen during the evolution of the system. 

Analysis of the algorithm in the infinite limit yields an exact solution for the flow 
in an infinite space. Simulated results agree perfectly with the analytic solution. 

Finite-space behavior is more complex because single sites are no longer indepen- 
dent. We have, however, been able to obtain for the average flow a non-trivial upper 
limit, which is exact ior m, k > L — 1. 
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